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Berezinskii-Kosterlitz-Thouless Transition of Spin-1 XXZ Chains 
in a Staggered Magnetic Field 
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Phase diagram of S = 1 XXZ chain in a staggered magnetic field is obtained numerically. The 
Berezinskii-Kosterlitz-Thouless transition between the XY and the antiferromagnetic phases is 
studied by the level-spectroscopy. Moreover we find that there is no distinction between the 
Haldane and the antiferromagnetic phases, since the Gaussian critical line does not appear, in 
contrast to the 5=1/2 case. It is expected that some hidden Z2 x Z2 symmetry breaking 
happens not only in the Haldane phase but also in the antiferromagnetic one. 

KEYWORDS: XXZ chain, staggered magnetic field, Berezinskii-Kosterlitz-Thouless Transition, level-spectroscopy, 
hidden Z2 X Z2 symmetry 



§1. Introduction 

The groundstatc properties of the spin-1 XXZ chain 
have_been studied extensively since Haldane's conjec- 
ture.EP The XXZ chain is described by the following 
Hamiltonian, 

L 

H XXZ = J2( S ? S f+i + S f S f+i + A ^? +1 ), (1.1) 

i=l 

where = S(S + l) n Fbr S = 1/2 case, the groundstate 
was obtained exactly™ and its properties are understood 
well.BcP The ferromagnetic phase is for A < — 1 and the 
Neel phase for A > 1. The XY phase for -1 < A < 1 
is characterized by the gapless excitation and the power- 
law decay of the spin correlation function. Before 1983, 
one had thought that the phase structure for S = 1/2 is 
valid forJiigher S, from the spin-wave picture. However, 
HaldanetP predicted that a novel phase (A c i < A < A c2 ) 
appears between the XY phase and the Neel phase for in- 
teger S, in which there is an excitation gap and the spin 
correlation function decays exponentially, in contrast to 
the half-odd integer S case. He also suggested that the 
Berezinskii-Kosterlitz-Thouless (BKT) transitiorJj'B' oc- 
curs at A c i, and that the transition at A C 2 belongs to 
the 2D Ising universality class. Mo'Qiadays his prediction 
is supported by many authors.Er&BEIl'lIS' A number of 
numerical studies with respect to S = 1 have shown that 
A cl = ~ 0.3 and A c2 = 1T8±0.02.BBBBB§> Al- 
though it was difficult to estimate A c i precisely duc_to 
the BKT transition, pcently it is settled that A ci = 0.E3 
Kennedy and Tasakili3> found that the S =1 XXZ chain 
contains a hidden Zi x Zi symmetry, and that the nature 
in the gapful phase, or Haldane phase, is explained by 
the complete breaking of the Zi x Zi symmetry. 

The hidden Zi x Zi symmetry of the groundstate can 
be controlled by the bond-alternation, i.e. replacing the 
nearest neighbor interaction in cq.(l.l) with 1 — SL^l 
In the isotropic case A — 1, Affleck and HaldaneOJ'l 



argued that the bond-alternating system is mapped onto 
the 0(3) nonlinear a model with the topological term 
whose angle is 9 = 2^5(1 — S), and that the excita- 
tion becomes massless when 6 /it is odd integer. There- 
fore there should exist 2S transition points ancL2.5 + 1 
massive phases between — 1 < S < 1. OshikawaEJ sug- 
gested that, (although his argument is not applicable to 
half-odd integer S) the Zi x Zi symmetry is essential 
for the successive dimerization transitions. The mas- 
sive phases are classified into two types, depending on 
whether the Zi x Zi symmetry is spontaneously bro- 
ken or not, and such two kinds of phases appear alter- 
nately. For S = the— HaJdane-dimer transition oc- 
curs at S c = 0.2598,Byk4k§ and its universality is the 
same class as the level k — 1 SU(2) Wess-Zumino-Witten 
modelrJ) The whole phase diagram of the S = 1 bond- 
alternating system with arbitrary A was obtained by Ki- 
tazawa, Nomura and OkamotoBS' They pointed out that 
a Gaussian critical line lies between the Haldane and 
the dimer phases, on the analogy of the Ashkin- Teller 
model with the Zi x Zi symmetry. Moreover they stated 
that the XY-Haldane phase boundary is just on A = 0. 
The XY-Haldane and the XY-dimer transitions are of 
the BKT type. The Haldane-Neel and the dimer-Neel 
transitions belong to the 2D Ising universality class. 

In this paper, as another example with the staggered 
interaction, we investigate the BKT transition of the 
XXZ chain in the staggered magnetic field, whose Hamil- 
tonian is given by, 



# = ffxxz + A^(-irSf 



(1.2) 



The staggered magnetic field induces the antiferromag- 
netic (AF) order, which is not the spontaneously sym- 
metry breaking. When A 3> 1, the AF phase character- 
ized by the singlet groundstate, AF long-range-order and 
gapful excitation should appear. Recently, Alcaraz and 
Malvezzl^P examined this model for S = 1/2. They 



2 



Masayoshi TSUKANO and Kiyohide NOMURA 



determined the XY-AF phase boundary by the phc- 
nomenological renormalization $£©up (PRG) method. c3 
However, Okamoto and NomuracS criticized their study, 
since the XY-AF phase transition is thought to be of the 
BKT type, and in that case the simple application of the 
PRG method may lead to the false conclusion. To over- 
come such a difficulty, one of the autherscj) proposed a 
new powerful method, " level-spectroscopy" , to estimate 
the BKT transition point precisely. This idea is based on 
the SU(2)/Z2 symmetry at the BKT transition point. In 
this paper, we apply the level-spectroscopy to the XY-AF 
transition in the staggered magnetic field, and we obtain 
the phase diagrams for S — 1/2 (in Fig. 1) and for S = 1 
(in Fig. 2). In the 5 = 1/2 case, there are two BKT lines 
and one Gaussian critical line (— 1/V2 < A < 1, A = 0). 
The Gaussian line_can be determined by the recent Ki- 
tazawa's method.EP The Neel phase (A > 1,A = 0), in 
which the groundstate is doublet, constructs a first-order 
transition line. For S = 1, on the other hand, although 
the BKT line exists, any Gaussian line does not appear. 
This is quite different from the bond-alternation case. 
The Haldane phase (0 < A < A C 2,A = 0) is connected 
continuously to the AF phase. The Neel phase is on the 
line A > A C 2, A = 0. The existence of a hidden Z 2 x Z 2 
symmetry is discussed in summary. 

§2. Method and Symmetry 

We consider the 2D sine-Gordon model, which is 



derived from the Hamiltonian (1.2) by means of the 
bosonization technique,© as an effective theory to de- 
scribe the BKT transition. The action is written by 



S = 



2-kK 



J drdx(V<j)) 2 + ^^ J drdxcosv^, (2.1) 

where a is a short distance cut-off. We also define the 
dual field 9(t, x) as 

^(T,x) = -iK^9(r,x), (2.2a) 
d d 

—4>(T,x)=iK—9(T,x). (2.2b) 

Although the U(l) symme try for the field 4> is broken 
by the second term of cq.( [2.l| ), it remains in the dual 
field 9. Here we compactify <f>, 9 on a circle with radius 
1/V2- The winding number for a configuration of 4> is 
nothing but total magnetization along the z-axis. The 
coupling constants K and y^ are renormalized by the 
scaling transformation a — > ae d ' n obeying the following 
renormalization group equations! 



dl 
d 
dl 



uo(i) = -$(0, 
vS) = -yo(i)y<t>{i), 



(2.3a) 
(2.3b) 



where yo = K/2 — 2. For the finite system, I is related 
to L by I = logL. The renormalization flow diagram 
is sketched in Fig. 3. The perturbation of cos\/20 is 
irrelevant for K > 4 and the trajectory flows into the 
Gaussian fixed line (yo > 0,y<p = 0), which corresponds 
to the XY phase with no symmetry breaking. When 
K < 4, cos V2<p is a relevant operator and y^ goes to 



infinity. In this case, ((f)) — tt/^/2 and the translational 
symmetry is broken. This means that the groundstate 
has AF long-range-order caused by the staggered mag- 
netic field. The BKT line (y^ = yo) is the boundary 
between the XY and the AF phases, and K is renormal- 
ized to 4. The scaling dimension of the vertex operator 
O n ,m — exp(m v / 20) exp(imv / 2^) on the Gaussian model 
(y<f> = 0) is given by 



n 2 K 



m 



(2.4) 



Thus the scaling dimension of cos V2(j) becomes 2 at the 
BKT fixed point (y^ = yo = 0). By the way, besides 
cos V2(j), there are four operators which become marginal 
at the BKT fixed point. They are sin V^cj), exp(±i4V26>) 
and the marginal operator defined by 



M = —{v<t>y. 



(2.5) 



However, this degeneracy splits by the logarithmic cor- 
rections on the BKT line away from the BKT fixed point 
(y<P = 2/o = 0). Next we consider how such splitting hap- 
pens. Here we note that M and cos ^/2<f) are hybridized 
near the BKT line, 



A = M + cos V2(j), 

B = V2cos\/20- -?-=M. 

V2 



(2.6) 
(2.7) 



The coefficients arc determined by the orthogonal con- 
dition (A(ti, Xi)B(t 2i x 2 )) = 0. The renormalized 
scaling dimensions of the operators, A, B, sin \f2<\> and 
exp(±i4y26>) are given as, up to the lowest order in y , 



(2.8a) 



ar (0=2- tt ,(0(l+gi), 



Xl (l)=2 + 2y (l) [l + -t 
x 2 (l)=2 + y (l), 



(2.8b) 
(2.8c) 



x 3 (l)=2-y {l), (2.8d) 

respectively.il' An additional parameter t plays the role 
of the deviation from the BKT line, y^/yo = 1 + t, and 
yo(l) = l^ 1 = (logL) -1 when (< 1. As an application 
of the conformal field theory (CFT), it is known that the 
scaling dimension-is related to the finite-size gap of the 
periodic system,EJ 

2-77" 7? 

E n (L) - E g {L) = — x n , (2.9) 

where v denotes the spin-wave velocity. Hence the eigen- 
values corresponding to xo and £3 cross linearly on t 
and degenerate on the BKT line, reflecting the SU(2)/Z2 
symmetry of the BKT line. 

We apply the above sine-Gordon theory to the orig- 
inal spin system. The Hamiltonian (1.2) is invariant 
under spin rotation around the z-axis, translation by 
two sites (Si — > Si +2 ) and space inversion (Si — > 
Sl-i'-, reflection on a spin site). Therefore eigenstates 
are classified by total spin moment (S^ = 
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wave number (q = 2im/L) and parity (P = ±1). Besides 
them, symmetry operation includes modified time rever- 
sal (T :' Sf -> Sf +1 ,Sf -> S? +1 ). The symmetry of 
the sine-Gordon operators is summarized in Table I. As 
a result, we can determine the BKT line by level-crossing 
of the excitation with q — 0, SJ, = 0, P = 1, T = 1 and 
that with q = 0, Sf. = 4, P = 1. 

Table I. Correspondence of the excitations to the sine-Gordon 
operators. 



Symmetry 


Operator on 


Notation of 


(9,5*,p,r) 


sine-Gordon model 


scaling dimension 


(0,0,+,+) 


M 


XO 


(0,0,+,+) 


cos \/24> 


X\ 


(0,0,-,-) 


sin \f2(j> 


XI 


(0,4,+, ) 


cxp(-i4V29) 


X3 



§3. Numerical Results 

The phase diagram of the S = 1 system is shown 
in Fig. 2. We demonstrate numerically the XY-AF 
transition, fixed A = —0.5 especially. Figure 4 shows 
some low-energy excitations of the finite-size system with 
L = 14. The critical point A c is determined by the cross- 
ing point of the excitation energy in q — 0, Sy. — 0, P — 
1, T = 1 subspace and that in q = 0, Sf. = 4, P = 1 sub- 
space. We obtain A C (Z) up to L = 18 (see Fig. 5), and 
then we estimate that A c = 0.855 ± 0.001 in the ther- 
modynamic limit L — > oo. The correction of 1/L 2 origi- 
nates from the irrelevant field with the scaling dimension 
4. At the critical point A = —0.5, A — 0.855 obtained 
above, we also calculate the averaged scaling dimensions 
defined as (xo + X2)/2 and (2xq + xi)/3, which elimi- 
nate the logarithmic correction, from the energy gaps in 
As is shown in Fig. 6, they converge to 



C qs.(y),(^9D 

2 within 1% error. The conformal anomaly is also esti- 
mated to be c = 0.993. These features confirm that the 
XY-AF transition is of the BKT type. The XY-AF tran- 
sition line approaches the point A = A = 0, and close to 
A = A = 0, it fits on A 2 oc —A very well (see inset in 
Fig. 2). 

To investigate the Gaussian critical line, we make use 
of the Kitazawa's mcthod.Ef According to him, when 
antiperiodic boundary condition (S^' v — —Sq' v ,S^ — 
Sq) is taken, two low- lying energies with = have 
to degenerate at a Gaussian critical point. For S = 1/2, 
we calculate the energies of the antiperiodic system with 
L = 20 and A = 0.5. As is shown in Fig. 7, they cross 
linearly at A — 0. This means that the Gaussian line 
is drawn-ppu A-.— and is governed by the c = 1 U(l) 
CFT.BByM As to S = 1, we perform the similar 
calculation about L = 14, A = 0.5. Unlike the S = 1/2 
case, the level-crossing is not observed in Fig. 8. Thus 
the groundstate does not undergo the phase transition 
between the Haldane and the AF phases. 

§4. Summary and Discussion 

We have studied the effect of the staggered magnetic 
field on the S = 1 XXZ chain. The BKT transition 



line between the XY and the AF phases has been ob- 
tained numerically by the level-spectroscopy. Next we 
have tried to find the Gaussian critical line, using the 
Kitazawa's method. However, the Gaussian line is ab- 
sent and there is no distinction between the Haldane and 
the AF phases. Then we suppose that the Hamiltonian 
has some hidden Z 2 x Z 2 symmetry as the sine- Gordon 
model does, and we expect that the Zi x Z 2 symmetry 
is spontaneously broken in both phases. In contrast, the 
5=1/2 case can be interpreted more naturally from the 
hidden Z 2 x Z2 symmetry picture; it is a variant of the 
Ashkin-Tellar model. One can consider that the hidden 
Zi x Z2 symmetry is fully broken for A < 0, unbroken 
for A > and partially broken (Neel region) on the line 
A = 0, A > 1. Note that the definition of the fully Z 2 x Z2 
broken region and the unbroken region is artificial, sim- 
ilar to the S — 1/2 bond-alternating system where one 
can interchange the ordered and the disordered region, 
taking the pair (2i — JL 2z) or (2i,2i + 1) in the string 
order parameter .EEIljE^ 

With respect to the S = 1 case, the existence of 
the hidden Zi x Zi symmetry shall be more clarified 
by taking account of uniaxial single-ion anisotropy, i.e. 
D Yli=i i^i ) 2 ' a P ar t from the question what is the ex- 
plicit form of the Zi x Zi symmetry. At least when 
D ^> 1 and A = 0, the Z2 x Z2 symmetry breaking does 
not occur because the groundstate is in the "large-!}" 
phase. EJ' Thus, if the hidden Z2 x Z 2 symmetry is con- 
tained in the Hamiltonian, the Gaussian critical line has 
to emerge between the AF and the large^D phases. We 
are now investigating these possibilities. c3 
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Figure Caption 



Fig. 1. Phase diagram of the S = 1/2 system. The XY-AF tran- 
sition is of the BKT type, and the dotted line denotes the Gaus- 
sian critical line. The Neel phase (doubly degenerate) makes the 
first-order transition line (bold line). 

Fig. 2. Phase diagram of the S=l system. The inset shows the 
blow-up region for —0.1 < A < on double logarithmic scale. 
We obtain that a = 0.52 ± 0.02 on assumption that A <x (— A) CT 
near A = A = 0. 



Fig. 3. Renormalization flow diagram of the 2D sine-Gordon 
model. 



Fig. 4. Low-energy excitations as L = 14 and A = —0.5. Eigen- 
states belong to q = 0, Sf, = 0, P = 1, T = l(o), q = 0, Sf, = 
0, P = -1,T = -l(o) and q = 0, fif, = 4, P = 1(+) subspace, 
respectively. 



Fig. 5. Estimation of the BKT transition point as A = —0.5. 
Using the system sizes L = 8 ~ 18, we obtain that A c = 0.855 ± 
0.001. 



Fig. 6. Scaling dimensions, (xq + X2)/2(o) and (2xq + xi)/3(o), 
at the critical point A = —0.5, A = 0.855. 

Fig. 7. Low- lying energies in S^, = subspace of the 5 = 1/2 
system with L = 20 and A = 0.5, on antiperiodic boundary 
condition. 



Fig. 8. Low-lying energies in = subspace of the S = 1 
system with L = 14 and A = 0.5, on antiperiodic boundary 
condition. 
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